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ABSTRACT 

We report the discovery of a near-infrared (nIR) counterpart to the persis- 
tent neutron-star low-mass X-ray binary 4U 1705—440, at a location consistent 
with its recently determined Chandra X-ray position. The nIR source is highly 
variable, with i^^-band magnitudes varying between 15.2 and 17.3 and addi- 
tional J- and if-band observations revealing color variations. A comparison with 
contemporaneous X-ray monitoring observations shows that the nIR brightness 
correlates well with X-ray flux and X-ray spectral state. We also find possible 
indications for a change in the slope of the nIR/X-ray flux relation between dif- 
ferent X-ray states. We discuss and test various proposed mechanisms for the 
nIR emission from neutron-star low-mass X-ray binaries and conclude that the 
nIR emission in 4U 1705—440 is most likely dominated by X-ray heating of the 
outer accretion disk and the secondary star. 

Subject headings: accretion, accretion disks — binaries: close — stars: individual 
(4U 1705-440) — X-rays: binaries 



^This paper includes data gathered with the 6.5 meter Magellan Baade Telescope located at Las Campanas 
Observatory, Chile, and the 4 meter Blanco Telescope located at CTIO, Chile. 
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Introduction 



4U 1705—440 is a neutron-star lo w-mass X-ray binary (NSXB) that was discovered in the 
early 1970's with the Uhuru satellite (iGiacconi et al.lll972l ). The X- ray spectral and variabil- 
ity pr operties place it among the class of the so-called atoll sources (IHasinger fc van der Klis 
19891 ). persistent NSXBs that accrete between a few and a few tens of perce nts of the Ed- 
dington luminosity. The s ource shows type-I (therr nonuclear) X-ray bursts (jSztajno et al. 
1985 : Gottwald et al. ^989), strong ra pid variability ( Olive et al. 2003 ). including kHz quasi 



periodic oscillations ( Ford et al. 19981). and some of the strongest long-term X-ray variations 
seen in persistent NSXBs ( iReig et al.ll2003l ) . These long-ter r n vari ations correspond to X-ray 
spectral-state transitions, as was shown by lBarret fc Olivd (120021 ). 



Despite the fact that 4U 1705—440 has been known for more than 35 years, no optical, 
infrared, or radio counterparts have been identified. Given the pronounced X-ray state 
transitions and large luminosity changes in 4U 1705—440, such counterparts are expected 
to be variable as well, which should in principle facilitate identification efforts. The main 
reason that counterparts have not been found is that, until recently, no accurate X-ray 
source position was available, impeding identification efforts in the c r owde d Galactic-bulge 
field. This changed with a Chandra observation by Idi Salvo et al.l (120051). which yielde d 
source coordinates with an absolute uncertainty of ~0'.'5. As noted by Idi Salvo et al.l (120051 ). 
these new coordinates differ significantly from previously reported ones, by ~0'15. 

We have conducted a search for the near-infrared (nIR) counterpart of 4U 1705—440, 
taking advantage of the improved Chandra position. The main benefit of searching in the 
nIR in compar ison to the optical i s the much lower extinction; based on Nh estimates from 
X-ray spectra (jPiraino et al.ll2007l) we expect A y ^ 10.8, but only Ak ~ 1.3. Moreover, as 



was recently shown by lRussell et al.l (120061 . 120071 ). the nIR can provide important constraints 



on the emission properties of the jet and accretion disk in NSXBs. In particular, for atoll 



sources iRussell et al.l (120071 ) found that above an X-ray luminosity of 10^^ ergs~^ the nIR is 
dominated by jet emission, while below these luminosities X-ray reprocessing in the accretion 
disk dominates (although thermal emission from the disk cannot be ruled out). 

In this paper we report the discovery of a variable nIR counterpart to 4U 1705—440. In 
^we present our nIR observations along with a set of quasi-simultaneous X-ray observations. 
The identification of the nIR counterpart and a study of nIR/X-ray correlations are presented 
in 1131 Finally, in §11 we discuss our results in the framework of several proposed mechanisms 
for nIR emission in NSXBs. 
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2. Observations &: Data Reduction 

2.1. near-Infrared 

We observed the field of 4U 1705—440 a number of times in the nIR between 2006 
May 16 and 2007 September 22. Our primary observations were with the 6.5 m Magellan I 



(Baade) telescope and Persson's Auxiliary Nasmyth Infrared Camera (PANIC; iMartini et al. 



20041 ). which contains a 1024 x 1024 HAWAII I HgCdTe detector with a 0'.'125 plate scale 
for a 128" x 128" field of view. We also obtained two observations with the 4 m Blanco 
Telescope at the Cerro Tololo Inter- American Observatory (CTIO), using the Infrared Side 
Port Imager (ISPI). ISPI has a 2048 x 2048 HgCdTe HAWAII-2 array with 0"3 plate scale, 
for a 10^25 x 10^25 field of view, although to match the PANIC data and to avoid dealing 
with geometric distortions we only used the central 512 x 512 pixels. The majority of the 
data were taken with a K^^ filter, but we did obtain some observations with the J and H 
filters. A log of the observations is given in Table [TJ There we list the total exposure times 
for each observation, although the data actually consist of a number of short (10-20 s for 
i^s; and up to 48 s for the other filters) exposures taken at positions dithered by up to 10". 

We reduced the data using custom routines in PyRAF. First we linearized the data using 
the coefficients given on the PANld and ISPfl web sites. We then subtracted dark frames, 
divided by a fiat field (using twilight fiats for PANIC and dome fiats for ISPI), and then 
produced a sky frame for subtraction by averaging the four exposures on either side of a 
given exposure. We then added the sky-subtracted exposures together, identified all the 
stars, and produced masks for the stars to improve the sky frames in a second round of sky 
subtraction. The coaddition of the individual exposures shifted the images of all the bands 
from a given observation to the same position. 

For a reference observation, we use d Epoch 1 from Table H] (2006 May 16). We deter- 



mined the astrometry using 64 2MASS (ISkrutskie et al.ll2006l ) stars that were well isolated 



and not very saturated, fitting for a plate scale, rotation, and zero point. The solution for 
the Ks image had an rms of 0'.'07 in each coordinate, and we transferred that solution to the 
other bands. For photometry, we used 34 2MASS stars that were unsaturated in all bands, 
well detected in 2MASS, and well isolated in all the data. We determined a zero point for 
each band (not including any color terms), and got zero-point uncertainties of 0.09 mag. 

All of the other data were corrected to the astrometry and photometry of Epoch 1. 



" http : //www. Ico . cl/lco/telescopes-inf ormation/magellcin/instruments-1/panic/panic-documentation/paiiic- 
' http : //www . ct io . noao . edu/ instruments/ ir jLnstruments/ ispi/New/UsersGuide/linearity . html 
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A list of star positions an d magnitudes for 820 st ars that were well-detected was produced 
by running sextractor (IBertin &: Arnoutd Il996l ) on the stacked images of each band of 
Epoch 1. We then ran sextractor to measure each of those stars in all of the other PANIC 
images. To compute the astrometric solution we used about 450 stars in each exposure 
and got rms residuals of ~ O'.'OS in each coordinate. We then transferred the solutions to the 
J and H bands. The photometric referencing was done using 16 stars within 40" and ±2 mag 
of 4U 1705—440 (based on the 2006 May 16 observation; see below), shown in Figure [H We 
selected these stars by searching through the 820 initially identified stars for those with the 
lowest variation between all of the exposures, after correcting for different zero points. These 
stars had rms variations of < 0.05 mag over all the observations. 



2.2. X-rays 



During part of the period in which the nIR observations of 4U 1705—440 were done, 
the source was monitored in X-rays with a number of instruments. For this paper we m ade 
use of publicly available data from the All-Sky Monitoilfl (ASM; Levine et al. 1996 ) on- 
board the Rossi X-ray Timing Explorer (RXTE) and the Swift Burst Alert Telescop^ (BAT; 
Barthelmy et al.ll2005l ). The ASM and BAT typically observed the source a few times per 
day in the 2-10 keV and 15-150 keV bands, respectively. In the period between 2007 Aug 11 
and 2007 Sep 11 no BAT data were available due to attitude control problems on the Swift 
spacecraft 

A dditional X-ray dat a from pointed observations with the Proportional Co unter Array 
fPCA: ljahoda et al.lEooeh and High Energy X-ray Timing Experiment (HEXTE; iRothschild et al 
19981 ) onboard RXTE were analyzed. These instruments provide coverage in the 2-60 keV 
and 15-250 keV bands, respectively. Starting from 2007 June 26 we observed 4U 1705—440 
on average once every four days. The PCA and HEXTE data were analyzed with FTOOLS 
V6.3. We obtained background-corrected count rates from the PCA and HEXTE in energy 
bands that closely match the energy bands of the ASM and BAT: 3-10 keV for the PCA 
and 15-150 keV for HEXTE. 



The PCA and HEXTE data were also used for a spectral analysis. For the data reduction 
and modeling o f the spectra we followed the data-reduction and analysis steps taken by 
Lin et al.l (120071 ). Spectral fits were made in the 3.5-150 keV band, with the fixed to 



'■^ ht t p : / /xte . mit . edu/ ASM_lc . html 

^http:/ /heasarc. gsfc.nasa.gov/docs/swift/results/transients/ 
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1.9 X 10^2 atoms cm-2 flPiraino et al.lboOTh . 



Results 



3.1. Near-Infrared Counterpart Identification 



Examining our Epoch 1 data, we find a single bright object within the Chandra error 
circle (J2000): a = 17''08'°54!47, 6 = -44°06'07M, with an uncertainty (90%) of 0'.'6 de- 
termined by the Chandra absolute pointing uncertainty (see Fig. [1]). This nIR source, at 
a = 17'^08™54!48, 6 = — 44°06'07'.'4, is well within the Chandra error circle, and is consistent 
with an unresolved point source in all of the data. In the Epoch 1 observations, this source 
has K, = 15.23 ±0.02, H = 16.03 ±0.02, and J = 17.14 ±0.02, where the errors only reflect 
the measurement uncertainties. 

When comparing all observations, it becomes apparent that the source varies consid- 
erably. These variations are not caused by the different seeing conditions. In Figure [2] we 
plot all our photometric data, with the photometry listed in Table [21 Uncertainties of ~0.05 
mag arising from our flat-fielding and sky-subtraction procedures were not included in the 
Table, but were taken into account in all figures. We include the two best comparison stars 
in Figure [2], which have rms brightness variations of < 0.06 mag. In the Kg-hand, in which 
we have the most observations, the counterpart varies by as much as 2 mag. Similar varia- 
tions are seen in the other bands, although there are considerable color variations; H — Kg 
varies by ~0.4 and J — i^'^ by ~0.6. These color changes appear to be uncorrected to 
changes in nIR brightness. In Figure [3] we plot a few representative nIR spectral energy 
distributions. The data for t he spectral energy d istributions were de-reddened assuming an 



Nh of 1.9x10^2 atoms cm ^ (jPiraino et al. 



fact ors at particu l ar wa velengths, following 
and IPutra eTaP J2002h : Ak, = 1.27, Ah 



20071). which was convert e d to various reddenin g 



Predehl fc SchmittI fll995h . lCardeUi et all (119891 ) . 
2.04, and Aj = 3.03. A comparison with the 
slopes of disk and jet spectra, which are also shown in Figure |3l seems to rule out significant 
optically thin jet emission in the nIR (see also Q . 



3.2. Near-Infrared/X-ray correlations 

The X-ray light curves from the ASM (one-day averages), the BAT (four-day averages), 
and PCA/HEXTE (single observation averages) are shown in Figure HI The times of the nIR 
observations are indicated in the panel with ASM data and the corresponding nIR epoch 
numbers are given on top. The X-ray state-transition cycles can clearly be seen, with the 
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15-150 keV count rates increasing during periods of low 2-10 keV flux. The nIR properties 
of 4U 1705—440 were studied as a function of X-ray flux and as a function of position in 
X-ray color diagrams. 



3.2.1. X-ray fluxes 

For the X-ray flux study we only made use of PCA and HEXTE data, since the ASM and 
BAT data are not suitable for spectral fitting. We only used PCA/HEXTE data that were 
taken within two days of our nIR observations; this resulted in 10 observations, two of which 
were corabined (for nIR Epoch 10). The PCA and HE XTE spectra wer e fit in XSPEC VI 1.3 
krnaudlEggeh with simple phenomenological models (ILin et al.l 120071 ). For spectrally-hard 



observations this model consisted of a black body, to represent the boundary layer between 
the disk and neutron star, and a broken power law, for the non-thermal emission. For 
spectrally-soft spectra an additional disk black body was needed, to fit the accretion disk 
component. The models also included an absorption component that was fixed at 1.9x10^^ 
atoms cm~^. These models resulted in good fits for all our observations. Unabsorbed 3-100 
keV fiuxes were extracted and in Figure [5] we plot the i^'s-band fiux density versus the 3-100 
keV X-ray fiux; error bars for the X-ray fiuxes were set at 5%. A clear, positive correlation 
can be seen; fits with a power law resulted in a nIR/X-ray slope of 0.66±0.02 (see dashed 
line in Fig. [S]). Separating the data based on relative strength of the thermal (blackbody and 
disk black body) and non-thermal (power-law) components in the X-ray spectrum, suggests 
that the nIR/X-ray fiux relation is considerably steeper for non-thermal dominated (3-100 
keV fiux contribution > 75%) X-ray spectra (1.33±0.11; open squares) than in the thermal 
dominated (> 85%) ones (0.37±0.03; filled circles), however, we note that the number of 
points in both selections is small and that the non-thermal data points span a very narrow 
range in X-ray fiux. We also point out that the low energy boundary of 3 keV means that a 
large fraction of the fiux falls outside our fiux band, especially for the soft spectra. 



3.2.2. X-ray color diagrams 

Spectral states of NSXBs are often studied in terms of color-color or color-intensity 
diagrams. Since the count rates from the ASM, BAT, and HEXTE were often too low to 
produce reliable colors from them, our diagrams were created from PCA data only. The PCA 
data were supplemented by publicly available data of 4U 1705—440 from RXTE observation 
cycles 9 and 10 (February 2004 until November 2006), allowing us to put our own observa- 
tions in a better context. The color-color and color-intensity diagrams are shown in Figure El 
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in which we also identi fy the spectrally-hard ( "extreme isla nd" ) , transitional ( "island" ) , and 
soft ("banana") states (ILin et al.l 120071 : Ivan der Klisll2006l ). These spectral states appear as 
well-separated branches in the color-color diagram. Since only nine of our nIR observations 
could be matched with PCA data we explored an alternative way to identify spectral states 
that could also be applied to the ASM, BAT, and HEXTE data of 4U 1705-440. As a 
solution we simply plotted count rates from a low energy band (~2-10 keV) versus count 
rates from a high energy band (~15-150 keV). This was done for the ASM and BAT four-day 
averaged count rates and separately for our PCA and HEXTE data. To allow comparison 
between the ASM/BAT and PCA/HEXTE count-rate diagrams, the count rates of the var- 
ious instruments were renormalized. First, the ASM and BAT count rates were scaled by 
their maximum four-day averaged values (from at least 3a detections) from the period after 
the start of the BAT coverage (25.6 counts s~^ and 2.1x10"^ counts s~^, respectively). The 
scaling factors for the PCA and HEXTE light curves were determined visually by finding 
the best overlap with the renormalized ASM and BAT light curves, resulting in scaling fac- 
tors of 869 counts s^^ for the PCA and 17.4 counts s"^ for HEXTE. The scaled ASM/BAT 
and PCA/HEXTE count rate diagrams are shown in Figure [71 Both the ASM/BAT and 
PCA/HEXTE diagrams reveal similar two branched structures. The symbols/colors for the 
PCA/HEXTE panel are the same as for Figure [6] and a comparison with that figure shows 
that the left-side branch corresponds to the spectrally-hard ( "extreme island" ) state and the 
right-side branch corresponds to the soft ("banana") state. The transitional ("island") state 
is not well separated in these diagrams, largely overlapping with the hard state branch. Note 
that the branches of the ASM/BAT track are slightly shorter than those from PCA/HEXTE, 
mainly because of missing Swift data in the August /September 2007 period. 

To understand how the broadband X-ray flux changed along the tracks in Figure [7] we 
fitted the spectra of all the PCA/HEXTE observations in that figure (109 in total) using 
the same model as above and obtained unabsorbed fluxes in the 0.5-100 keV band, by 
extrapolating well below our low-energy boundary. The main reason for extrapolating down 
to 0.5 keV was to capture more flux from the thermal components, although uncertainties are 
involved in extrapolating to such low energies. We found that on the spectrally-hard branch 
the 0.5-100 keV flux correlated well with the HEXTE count rate and that on the spectrally 
soft branch the 0.5-100 keV flux correlated well with the PCA count rate. This suggests that 
mass-accretion rate increases as the source moves up either branch. The 0.5-100 keV flux 
range observed on the hard branch was 1.1-8.5x10"^ ergs~^cm~^; on the soft branch we 
measured 1.8-24x10"^ ergs~^ cm~^. Considering only the ten X-ray observations that were 
associated with the nIR observations, we find that the ones on the spectrally-soft branch all 
had higher 0.5-100 keV fluxes that the ones on the spectrally-hard branch. 

Using the same PCA observations as for our spectral analysis ( §3.2. ip . we overplot the 
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nIR data (grey circles) in the color-color and color-intensity diagrams in Figure O with the 
size of the symbols scaling linearly with the nIR flux density. Clearly, the brightest nIR ob- 
servations occurred in the spectrally-soft state, in which the X-ray flux reaches its maximum. 
We also selected ASM/BAT data that were closest in time to our nIR observations, using 
three- or four-day ASM/BAT averages that were centered on the times nIR observations. 
These data are plotted in Figure [TJ together with the PCA/HEXTE data (for which we used 
the same nIR/X-ray combinations as for Figures [5] and [6]) . Clearly, a large part of the full 
X-ray spectral range was covered by our nIR observations. The locations of the nIR epochs 
that were covered by ASM/BAT as well as PCA/HEXTE observations were consistent be- 
tween the two diagrams, except for Epoch 18, for which we derived a scaled ASM count rate 
of 0.96(3) counts s~^, but a scaled PCA count rate of 1.13(1) counts s~^. As can be seen 
from Figure [71 the additional ASM/BAT data strengthen our initial finding from Figure E] 
that nIR fluxes are highest in the spectrally-soft state. In Table [2] we have listed the spectral 
branch on which each of the nIR observations occurred; two observations occurred near the 
vertex of the two branches and identification is therefore uncertain. Comparison of Figure [7] 
and Table [2] shows that all soft-state observations have higher nIR fluxes than the hard-state 
observations. Since this is also true for the X-ray fluxes, the spectral-state dependence of 
the nIR flux might therefore simply reflect the nIR dependence on X-ray flux (see also §11). 



4. Discussion and conclusions 

We have made nIR observations of the field of the NSXB 4U 1705—440 and discovered 
a variable counterpart at the location of the Chandra X-ray source. The i^^-band magni- 
tude varied between 15.2 and 17.3 and additional J- and if-band observations revealed color 
changes that did not correlate with nIR brightness variations. Using 18 epochs of nIR obser- 
vations in combination with X-ray monitoring observations, we were able to study, for the 
first time, the nIR changes along an atoll color-color track. We found that 4U 1705—440 has 
a significantly higher nIR flux in the spectrally-soft ( "banana" ) state than in the spectrally- 
hard ( "extreme-island" ) state, although this might simply be the result of the higher X-ray 
fluxes (and mass-accretion rates) in the former state. The nIR flux density correlated well 
with X-ray flux, with possible indications for a steeper dependence in the spectrally-hard 
state. 

Comparing the monochromatic nIR luminosity of 4U 1705—440 (log(LxJ ~ 19.2 — 20.0 
[7.2 kpc], ~ 19.5 — 20.3 [9.6 kpc] ; for distance estimates, see below) with that of other NSXBs 



published in lRussell et al.l (120071 ). we find values that are very similar to syste ms with com- 
parable X-ray luminosities. For -Fs-ioofcevZ-^i^s we find values of 1.4-4.4x10^. iRussell et al. 
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( 120061 . 120071 ) discussed possible mechanisms for the nIR emission in NSXBs: optically-thin 
emission from the inner regions of a jet outflow close to the neutron star, thermal emission 
from an X-ray or viscously-heated disk, or thermal emissio n from the corapanion star. X-ray 
heating of the secondary was not explicitly considered by iRussell et al. J200d. l2007h 



even 



though it can be a significant source of optical emission, as shown by IZurita et al.l (120001 ). 
and presumably also of nIR emission. However, since it is difficult to distinguish heating of 
the disk from heating of the secondary with our data, we will not discuss the latter option 
separately. 

We start with an assessment of the contribution from a non-heated secondary star. The 
mass donors in low-mass X-ray binaries are typically late-type dwarfs or giants. We com- 
pare the absolute Kg magnitude s of 4U 1705—440 wit h published K values for K4 and M4 



dwarfs, and K5 and M 5 giants (IBessell fc BrettI Il988l : iBesselll Il99ll : iMikami fc Heckl Il982 



Carroll fc Os"tli3 Il996h . calculated using the distance estimates derived at the end of this 
section and a correction for absorption that is based on an Nh of 1.9x10^^ atoms cm~^. 
For dwarfs we find that they are three to seven magnitudes fainter than our lowest ob- 
served values. The absolute magnitudes of the giants, on the other hand, are three to 
four magnitudes brighter than our highest observed values. Note that these comparisons 
suffer from considerable distance and reddening related uncertainties. Nevertheless, these 
findings, in combination with the amplitude of the nIR variations, appear to rule out that 
thermal emission from a non-heated late-type dwarf or giant dominates the nIR emission of 
4U 1705—440. Given that the expected nIR emission from giants is higher than what we 
observe, the secondary in 4U 1705—440 is most likely a dwarf. This would then imply an 
orbital period of ~1-10 hrs for a Roche-lobe filling primary, which needs to be tested with 
follow-up photometry. 



T o disti nguish between the other possible mechanisms that are listed above, iRussell et al. 



(120061 . 120071 ) employed nIR/X-ray luminosity correlations and spectral energy distributions. 
In Figure [3] we compare the nIR spectral energy distributions of 4U 1705—44 with rnodels 
for optically-thin jet synchrotron emission (oc average slope taken from iFenderl 12006 



solid grey line), t he Ray leigh- Jeans tail of a multi-color disk black-body spectrum (oc z/^; 
Frank et al.lll992l . dashed line) and the v^/^ slope of a multi-color disk spectrum that is char- 
acteristic for frequencies higher than the Rayleigh- Jeans part (dotted line). The frequency 
of the break between the and u^f^ slopes of a multi-color disk black-body increases with 
mass-accretion rate and decreases with the size of the accretion disk. Slopes for an X-ray- 
heated accretion disk are not included in Fig ure [21 Such heating results in an additional 
bump in the optical/UV (jVrtilek et al.l Il990l ) and also extends the high-frequency end of 
the Rayleigh- Jeans tail to higher frequencies. Based on the model slopes in Figure [3l a 
significant contribution to the nIR from optically-thin jet synchrotron emission can be ruled 
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out; the slope of the nIR spectrum is the opposite of that expected from a jet and the 
relatively small changes in the shape of the spectrum also rule out significant changes in 
the nIR contribution from a jet across the observed nIR lu minosity range. Recen t Spitzer 
observations of another atoll source NSXB, 4U 0614+091, bv lMigliari et ali J2006h . suggest 
that one has to observe at lower IR frequencies to detect clear spectral signatures from a 
jet in NSXBs. In that source the jet synchrotron emission had a slope of —0.57, indicative 
of optically-thin emission and close to the —0.6 slope we adopted in Figure [31 The depen- 
dence of nIR flux on X-ray spectral state also argues against a jet interpretation of th e nIR 
flux. Observations of black- hole LMXBs (IHoman et al.l l2005l : iBuxton fc Bailynl |2004J ) have 
shown that the nIR contribution from a jet changes significantly during transitions between 
spectrally-hard and soft X-ray states, with the nIR contribution in hard states being higher. 
In 4U 1705—440 we observe the opposite behavior, suggesting that the nIR flux is dominated 
by other components. 

The slopes of the spectral energy distributions in Figure [3] are consistent with those of 
multi-color disk black body models, being closer the Rayleigh- Jeans tail on the low-frequency 
end and closer to the flatter part of a disk spectrum on the high-frequency end. We note that 
in recent literature on 4U 1705—440 Nh values in the range of 1.2-2.4x10^^ atoms cm~^ can 
be found (we adopted 1.9x10^^ atomscm"^). The low value of this range brings the slope 
closer to that of the flat part of a disk spectrum, but still far from the jet slope, whereas 
the high value steepens the data to bring them closer to the slope of the Rayleigh- Jeans tail. 
Based on the spectral energy distributions we are not able to distinguish between X-ray- 
heated or viscously-heated disk spectra. NIR/X -ray flux correlations could in principle be 
used to shed more light this. iRussell et al.l (120071 ) used such correlations for a large ensemble 
of NSXBs that covered more than five decades in X-ray luminosity. Although our X-ray 
observations of 4U 1705—440 only span about a decade in X-ray luminosity, a clear correlation 
between X-ray and nIR fluxes can be seen in Figure O The best-fit power-law slope of 
0.66±0.02 is consistent with that found for the large sample of NSXB, which d id not include 



4U 1 705—440, and also consistent with the expected slope from X-ray h eating (Russell et a 



20071 ). A viscously-heated disk would have resulted in a slope of ~0.3. iRussell et al.l (120071 ) 
combined the NSXB data regardless of spectral state. Dividing our data based on the 
dominant X-ray spectral component (thermal vs. non-thermal) reveals a significantly steeper 
relation for the (hard) non-thermal X-ray state than for the (soft) thermal X-ray state (see 
Figure [5]): 1.33±0.11 compared to 0.37±0.03. Although these different slopes are obtained 
from small numbers of points, they are both no longer consistent with X-ray heating. The 
slope for the (soft) thermal state is marginally consistent with viscous heating, while the 
slope for the (hard) non-therm al state is similar to that found between radio flux and X-ray 



flux in NSXBs, Lradio oc Lx ' (IMigliari &: Fenderl 120061 ) . For this slope also to be observed 
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in the nIR, the nIR-radio spectral index would have to be very fiat, requiring the nIR to 
originate from the optically-thick part of a jet synchrotron spectrum, for which we see no 
evidence however. 

An additional test to distinguish between X-ray and viscous heating is to search for 
lags between changes in the nIR and X-ray emission of a few days or more. Such lags 
are not expected when X-ray heating is important, as reprocessing of the X-rays is near- 
instantaneous. If, on the other hand. X-ray heating is not important, variations in the 
X-rays from the inner disk will lag variations in the nIR emission from the non-irradiated 
outer disk, since changes in the mass-accretion rate travel through th e disk on the viscou s 



time scale. Such delays have been observed in GX 339-4 for example (iHoman et al.ll2005l ). 
where a lag time of 15-20 days was measured. In Figure [H] we plot the broadband X-ray flux 
of the thermal component (disk and boundary layer; filled dots) , the non-thermal power-law 
component (open squares), and the Ks fiux density light curves from an X-ray fiux cycle 
in 2007 July-September that was sampled well in the nIR. The result is ambiguous, mainly 
because the source underwent a spectral transition in X-rays, making it hard to follow the 
disk component over a long period of time; changes in the disk X-ray fiux during a transition 
probably do not only reflect changes in the mass-accretion rate through the inner disk. 
However, the figure does show that the nIR light curve resembles that of the total X-ray fiux 
(thermal and non-thermal combined), suggesting the nIR fiux is produced by X-ray heating. 
Better sampled light curves would be needed to perform a more detailed study of time lags. 



Based on observed spectral indices iRussell et al.l (120071 ) concluded that in atoll NSXBs 



with 2-10 keV luminosities above 10^^ ergs ^ jet emission dominates the n I R- Distance es- 



timat es from radius expansion bursts (7.2-9.6 kpc lJonker fc Nelemansll2004l : iGalloway et al. 



20061 ). indicate that the 2-10 keV luminosity of 4U 1705—440 is well above this value. How- 
ever, our spectral indices are more consistent with disk emission. The observed 0.5-100 keV 
fiux range of 4U 1705—440 leads to a total luminosity range of 0.7 — 14.9 x 10^'' ergs~^ (7.2 
kpc) or 1.2 — 26.5 x 10^^ (9.6 kpc). These values are quite high for atoll sources and come 
close to those of the Z sources, possibly indicating a larger accretion disk than is typical for 
atoll sources. This could then explain why X-ray heating dominates the nIR in 4U 1705—440. 

We conclude that X-ray heating of the outer disk and secondary star is the most likely 
source of the nIR emission in 4U 1705—440. Detection of possible orbital modulations and 
obtaining a nIR spectrum might help in distinguishing the relative importance of these two 
mechanisms. More observations are needed to study a possible X-ray state dependence of 
the nIR/X-ray relation. 
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Table 1. Near-IR Observation Summary for 4U 1705-440 



Epoch Datet H J 







Time 


Exp. (s) 


Airmass 


Seeing 


Time 


Exp. (s) 


Airmass 


Seeing 


Time 


Exp. (s) 


Airmass 


Seeing 


1 


2006-May-16 


06:06 


50 


1.04 


0.6 


06:22 


240 


1.04 


0.7 


06:11 


450 


1.04 


0.7 


2 


2006-Jul-15 


00:34 


230 


1.12 


0.7 


00:47 


1296 


1.10 


0.7 


01:19 


1620 


1.07 


0.7 


3 


2006-Aug-04 


00:45 


180 


1.04 


0.6 


00:51 


864 


1.04 


0.7 


01:18 


1620 


1.04 


0.7 


4 


2006-Aug-08a 


00:45 


90 


1.04 


0.9 


00:07 


432 


1.05 


1.1 


00:20 


810 


1.04 


0.9 


5 


2006-Aug-08b 


23:47 


180 


1.06 


1.0 


23:54 


432 


1.06 


1.0 


00:06 


1215 


1.05 


1.0 


6 


2006-Aug-09 


23:37 


90 


1.07 


0.7 


23:43 


432 


1.06 


0.8 


23:55 


810 


1.05 


0.9 


7 


2006-Aug-26 


23:30 


270 


1.04 


0.7 


23:53 


1296 


1.04 


0.7 


00:35 


2430 


1.06 


0.8 


8 


2007-Apr-07 


09:42 


100 


1.05 


0.6 










09:35 


100 


1.04 


0.6 


9 


2007-May-09 


07:56 


300 


1.06 


0.7 


















10 . . . 


. 2007-Jul-05 


01:08 


200 


1.16 


0.8 










01:31 


200 


1.12 


0.7 


11* . 


. 2007-Jul-24 


01:53 


320 


1.03 


1.0 


















12* . 


. 2007-Jul-27 


06:31 


480 


1.96 


0.8 


















13 ... 


. 2007-Jul-31a 


01:45 


90 


1.04 


0.7 


01:36 


270 


1.04 


0.7 


01:27 


270 


1.04 


0.7 


14... 


. 2007-Jul-31b 


23:25 


110 


1.16 


0.7 










23:34 


270 


1.14 


0.7 


15 . . . 


. 2007-Aug-03 


02:06 


300 


1.04 


0.6 


















16... 


. 2007-Aug-20 


00:51 


180 


1.04 


0.6 










01:01 


360 


1.04 


0.5 


17... 


. 2007-Sep-Ol 


02:46 


220 


1.30 


0.4 


















18 ... 


. 2007-Sep-22 


01:38 


360 


1.35 


0.8 


01:29 


450 


1.32 


0.7 


01:23 


180 


1.30 


0.9 



tAU dates and times are UT. 

*A11 observations were done with Baade/PANIC, except for these which used Blanco/ISPI. 



-16- 



Table 2. Photometry of 4U 1705-440 





Ld-y 






A A fi crn i"i"n H 












u 

11 


7 

J 


1 


Soft 


15 


23 ±0.02 


16.03 ±0.02 


17.14 ±0.02 


2 


?* 


16 


18 ±0.02 


16.82 ±0.02 


17.96 ±0.02 


3 


?* 


16 


10 ±0.02 


16.93 ±0.02 


18.16 ±0.02 


4 


Hard 


16 


98 ±0.08 


17.41 ±0.03 


18.47 ±0.03 


5 


Hard 


16 


74 ± 0.03 


17.37 ±0.03 


18.39 ±0.03 


6 


Hard 


16 


83 ± 0.03 


17.42 ± 0.03 


18.48 ±0.03 


7 


Hard 


16 


66 ± 0.02 


17.34 ± 0.02 


18.39 ±0.02 


8 


Soft 


15 


51 ±0.02 




17.05 ± 0.02 


9 


Soft 


15 


77 ±0.02 






10 .... 


Soft 


16 


19 ±0.02 




17.94 ±0.03 


11 .... 


Hard 


17 


25 ±0.07 






12 .... 


Hard 


17 


10 ±0.06 






13 .... 


Hard 


16 


41 ± 0.03 


16.94 ± 0.02 


18.01 ± 0.02 


14 ... . 


Hard 


16 


54 ± 0.03 




18.23 ±0.03 


15 .... 


Hard 


16 


31 ± 0.02 






16 .... 


Soft 


15 


61 ± 0.02 




17.10 ±0.02 


17 .... 


Soft 


15 


37 ±0.02 






18 .... 


Soft 


15 


60 ± 0.02 


16.18 ±0.02 


17.14 ±0.02 



*The state identification of these two observations is uncer- 
tain. 

Note. — The quoted errors only reflect s extract or uncer- 
tainties. Additional uncertainties related to flat fielding and 
sky subtraction are estimated to be 0.05 mag. 
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Fig. 1.— Images of the field around 4U 1705-440 from the 2006 May 16 PANIC observa- 
tion. The main panel shows a 1'.5 x 1^5 field of the Kg image. The Chandra position of 
4U 1705—440 is marked with a 0'/6-radius circle. We also indicate with tick marks a number 
of the photometric reference stars (some are hidden behind the sub-panels or are off this 
image): the two most stable stars have crosses, while the rest have only ticks to the north 
and east. The sub-panels show 15" x 15" fields of the J (top left) and H (bottom left) 
images. North is up, and east to the left. 
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Fig. 2. — Photometry of 4U 1705—440. In the top panel we show the Kg (diamonds), H 
(squares), and J (circles) photometry of 4U 1705—440 plotted as a function of epoch number. 
In the bottom panel we show the variation in the magnitudes of the two most stable reference 
stars (see Fig. [1]). The colors/gray-scales represent the same filters, but here the two symbols 
(points and stars) correspond to the two reference stars. The date of each epoch is listed 
in the bottom panel. [(See the electronic edition of the Journal for a color version of this 
figure.] 
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Fig. 3. — Six representative nIR spectral energy distributions (SEDs) of 4U 1705—440, de- 
reddened assuming an Nh of 1.9x10^^ atoms cm~^ [Ay — 10.3). The nIR epoch of each 
of the SEDs is shown on the left-side of the Kg data points. The solid grey hne depicts a 
typical optically thin synchrotron spectrum from a jet (oc The dashed line represents 

the Rayleigh- Jeans tail of the spectrum of a viscously heated accretion disk (oc z/^) and 
the dotted line the shallower 'stretched' part that is characteristic of such a disk spectrum 

(oc I/V3). 



-20- 



near— IR Epoch 

3-6 11-15 17 

12 7 8 9 10 16 18 




100 200 300 400 500 
Days since 2006 May 1 



Fig. 4. — X-ray light curves of 4U 1705—440 from various instruments. Prom top to bot- 
tom: RXTE/ ASM, Swift/BAT, RXTE/PCA, and RXTE/BEXTE. The ASM and BAT hght 
curves show one- and four-day averages, respectively, whereas the PCA and HEXTE light 
curves were rebinned to one data point per observation. The times of the nIR observations 
are indicated by the vertical lines and numbers in and above the top panel. 
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Fig. 5. — Kg flux density versus unabsorbed X-ray flux (3-100 keV). The dashed hne shows 
the best power-law fit to the data with an index of 0.66(2). Filled circles correspond to 
observations in which the thermal components dominates the X-ray spectrum (> 85%), 
open squares ones in which the non-thermal flux dominates (> 75%). 
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Fig. 6. — Color-color (left) and color-count rate diagrams (right) of 4U 1705—440, created 
from RXTE/PCA data. Soft color and hard color were defined as count rate ratios in the 
bands (4.0-7.3 keV)/(2.3-4.0 keV) and (9.8-18.2 keV)/(7.3-9.8 keV), respectively. Typical 
error bars are shown on the right hand side of each panel. Triangles, squares, and open circles 
indicate the spectrally hard ("extreme island"), transitional ("island"), and soft ("banana") 
states, respectively. The gray circles depict observations coincident (within two days) with 
our nIR observations, with the size of the circles scaling linearly with nIR flux density. [(See 
the electronic edition of the Journal for a color version of this figure.] 
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Fig. 7. — High-energy versus low-energy scaled count rates from ASM/BAT (left) and 
PCA/HEXTE (right). The count rates were scaled as described in §3.21 Typical error 
bars are shown in the lower right corner of each panel. Triangles, squares, and open circles 
in the right panel indicate the spectrally hard ("extreme island"), transitional ("island"), 
and soft ("banana") states, respectively. The gray circles depict observations coincident 
(within two days) with our nIR observations, with the size of the circles scaling linearly with 
nIR flux density. [(See the electronic edition of the Journal for a color version of this figure] 
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Fig. 8. — Comparison of thermal (dots) and non-thermal (squares) X-ray light curves (3-100 
keV) and Kg-band (stars) light curves during a transition cycle in July-September 2007. 



